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1 Introduction 

In [Elj we introduced several computational challenges concerning Shimura 
curves, and some techniques to partly address them. The challenges are: ob- 
tain explicit equations for Shimura curves and natural maps between them; 
determine a Schwarzian equation on each curve (a.k.a. Picard-Fuchs equation, 
a linear second-order differential equation with a basis of solutions whose ratio 
inverts the quotient map from the upper half-plane to the curve) ; and locate CM 
(complex multiplication) points on the curves. We identified some curves, maps, 
and Schwarzian equations using the maps' ramification behavior; located some 
CM points as images of fixed points of involutions; and conjecturally computed 
others by numerically solving the Schwarzian equations. 

But these approaches are limited in several ways: we must start with a 
Shimura curve with very few elliptic points (not many more than the mini- 
mum of three) ; maps of high degree are hard to recover from their ramification 
behavior, limiting the range of provable CM coordinates; and these methods 
give no access to the abelian varieties with quaternionic multiplication (QM) 
parametrized by Shimura curves. Other approaches somewhat extend the range 
where our challenges can be met. Detailed theoretical knowledge of the arith- 
metic of Shimura curves makes it possible to identify some such curves of genus 
at most 2 far beyond the range of [El] (see e.g. [RoblGRj ), though not their 
Schwarzian equations or CM points. Roberts [Rob] showed in principle how to 
find CM coordinates using product formulas analogous to those of |GZ| for differ- 
ences between CM j-invariants, but such formulas have yet to be used to verify 
and extend the tables of |E1| . Errthum |Er| recently used Borcherds products to 
verify all the conjectural rational coordinates for CM points tabulated in |E1| for 
the curves associated to the quaternion algebras over Q ramified at {2, 3} and 
{2, 5}; it is not yet clear how readily this technique might extend to more com- 
plicated Shimura curves. The p-adic numerical techniques of |E3j give access to 
further maps and CM points. Finally, in the {2,3} and {2,5} cases Hashimoto 
and Murabayashi had already parametrized the relevant QM abelian surfaces 
in 1995 HMJ, but apparently such computations have not been pushed further 
since then. 

In this paper we introduce a new approach, which exploits the fact that some 
Shimura curves also parametrize K3 surfaces of Neron-Severi rank at least 19. 
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"Singular" K3 surfaces, those whose Neron-Severi rank attains the characteristic- 
zero maximum of 20, then correspond to CM points on the curve. We first en- 
countered such parametrizations while searching for elliptic K3 surfaces of max- 
imal Mordell-Weil rank over Q(i) (see [E4 ), for which we used the K3 surface 
corresponding to a rational non-CM point on the Shimura curve X(6, 79)/(i06-79) 
of genus 2. The feasibility of this computation suggested that such parametriza- 
tions might be used systematically in Shimura curve computations. 

This approach is limited to Shimura curves associated to quaternion algebras 
over Q. Within that important special case, though, we can compute curves 
and CM points that were previously far beyond reach. The periods of the K3 
surfaces should also allow the computation of Schwarzian equations as in [LY , 
though we have not attempted this yet. We do, however, find the corresponding 
QM surfaces using Kumar's recent formulas |Kuj that make explicit Dolgachev's 
correspondence |Doj between Jacobians of genus-2 curves and certain K3 surfaces 
of rank at least 17. The parametrizations do get harder as the level of the Shimura 
curve grows, but it is still much easier to parametrize the K3 surfaces than to 
work directly with the QM abelian varieties — apparently because the level, 
reflected in the discriminant of the Ncron Severi group, is spread over 19 Ncron 
Severi generators rather than the handful of generators of the endomorphism 
ring0 In this paper we illustrate this with several examples of such computations 
for the curves X(A, 1) and their quotients. As the example of X(6, 79)/ (1^6-79) 
shows, the technique also applies to Shimura curves not covered by X(A, 1), but 
already for X(A, 1) there is so much new data that we can only offer a small 
sample here: the full set of results can be made available online but is much too 
large for conventional publication. Since we shall not work with X(A, M) for 
M > 1, we abbreviate the usual notation X(A, 1) to X(A) here. 

The rest of this paper is organized as follows. In the next section, we re- 
view the necessary background, drawn mostly from |VilR ot2 BHPV], concern- 
ing Shimura curves, the abelian and K3 surfaces that they parametrize, and 
the structure of elliptic K3 surfaces in characteristic zero; then give A. Ku- 
mar's explicit formulas for Dolgachev's correspondence, which we use to recover 
Clebsch-Igusa coordinates for QM Jacobians from our K3 parametrizations; and 
finally describe some of our techniques for computing such parametrizations. In 
the remaining sections we illustrate these techniques in the four cases A = 6, 
N = 14, A = 57, and N = 206. For A — 6 we find explicit elliptic models for 
our family of K3 surfaces S parametrized by X(6)/ (wq), locate a few CM points 
to find the double cover X(6)— >X(6)/ (wq), transform S to find an elliptic model 
with essential lattice A css D £7 © E$ to which we can apply Kumar's formulas, 
and verify that our results are consistent with previous computations of CM 
points [Elj and Clebsch-Igusa coor dinates [HMj . For A = 14 we exhibit S and 
verify the location of a CM point that we computed numerically in [Elj but 
could not prove using the techniques of [E1IE3] . For A = 57, the first case for 

1 It would be interesting to quantify the computational complexity of such computa- 
tions in terms of the level and the CM discriminant; we have not attempted such an 
analysis. 
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which ~X.(N)/(wn) has positive genus, we exhibit the K3 surfaces parametrized 
by this curve, and locate all its rational CM points. For N — 206, the last case 
for which X(iV) / (wn) has genus zero, we exhibit the corresponding family of K3 
surfaces and the hyperelliptic curves X(206) and X(206)/(w2), X(206)/(wio3) 
covering the rational curves X(206)/(w206) and X(206)/(u>2, Wios)- 

2 Definitions and techniques 

Quaternion algebras over Q, Shimura curves, and QM abelian sur- 
faces. Fix a squarefree integer N > with an even number of prime factors. 
There is then a unique indefinite quaternion algebra A/Q whose finite ramified 
primes are precisely the factors of N. Let O be a maximal order in A. Since A 
is indefinite, all maximal orders are conjugate in A, and conjugate orders will 
be equivalent for our purposes. Let O* be the group of units of reduced norm 1 
in O; let T be the arithmetic subgroup OJ/{±l} of A*/Q*; and let r* be the 
normalizer of r in the positive-norm subgroup of A*/Q*. If N = 1 then r* = r; 
otherwise r* /T is an abelian group of exponent 2, and for each factor d\N there 
is a unique element Wd S -T* / T whose lifts to A* have reduced norms in d ■ Q* 2 . 

Because A is indefinite, A®qR is isomorphic with the matrix algebra M2CR), 
so the positive-norm subgroup of A*/Q* is contained in PSL 2 (R) and acts on the 
upper half-plane H.. The quotient Tt/r is then a complex model of the Shimura 
curve associated to r, usually called X(JV, f). In |Elj we called this curve X(\) 
in analogy with the classical modular curve X(l) (see below), since N was fixed 
and we studied Shimura curves that we called Xo(j>), Xi(p), etc., associated 
with various congruence subgroups of A*/Q*. In this paper we restrict attention 
to Ti/r and its quotients by subgroups of r*/T; thus we return to the usual 
notation, but simplify it to X(AT) because we do not need X(7V, M) for M > 1. If 
N = 1 then A = M 2 (Q), and we may take O = M 2 (Z), when r = T* = PSL 2 (Z) 
and Ti. must be extended by its rational cusps before we can identify TL/ T with 
X(l). Here we study curves X(AT) and their quotients only for N > 1, and these 
curves have no cusps. 

The Shimura curve X(iV) associated to a quaternion algebra over Q has 
a reasonably simple moduli description. Fix a positive anti-involution g of A 
of the form g(/3) = for some fx e O with fi 2 + N = 0. Then X(N) 

parametrizes pairs (A, l) where A is a principally polarized abelian surface and 
1 is an embedding of O into the ring End(A) of endomorphisms of A, such that 
the Rosati involution is given by g. See [Rot2l §2 and Prop. 4.1]. This gives X(iV) 
the structure of an algebraic curve over Q. 

An abelian surface with an action of a (not necessarily maximal) order in a 
quaternion algebra is said to have "quaternionic multiplication" (QM). A com- 
plex multiplication (CM) point of X(iV) is a point, necessarily defined over Q, 
for which A has complex multiplication, i.e. is isogenous with the square of a 
CM elliptic curve. We shall use the QM abelian surfaces A to find models for 
the Shimura curves X(AT) and locate some of their CM points. 
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When N — 1, an abelian surface together with an action of O = M2(Z) is just 
the square of an elliptic curve, so we recover the classical modular curve X(l). We 
henceforth fix N > 1. Then the group r* /T, acting on X(iV) by involutions that 
we also call Wd, is nontrivial. These involutions are again defined over Q, taking 
(A, i) to (Ad, id) for some Ad isogenous with A. Specifically, Ad is the quotient 
of A by the subgroup of the ci-torsion group A[d] annihilated by the two-sided 
ideal of O consisting of elements whose norm is divisible by d, and the principal 
polarization on Ad is 1/d times the pull-back of the principal polarization on A. 
In particular Ad is CM if and only if A is. Hence the notion of a CM point makes 
sense on the quotient of X(N) by r*/T or by any subgroup of r*/T. If a CM 
point of discriminant — D on X(N)/(r*/T) is rational then the class group of 
Q(y/—D) must be generated by the classes of primes lying over factors p\D that 
also divide N. Thus the class group has exponent 1 or 2 and bounded size; in 
particular, only finitely many D can arise. In each of the cases N = 6, 14, 57, 
and 206 that we treat in this paper, N has two prime factors, so the class number 
is at most 4 and we can cite Arno [Ar. to prove that a list of discriminants of 
rational CM points is complete. When N has 4 or 6 prime factors we can use 
Watkins' solution of the class number problem up to 100 |Waj . 

We have An = A as principally polarized abelian surfaces, but for N > 1 the 
embeddings i, l n are not equivalent for generic QM surfaces A. When we pass 
from A to its Kummcr surface we shall lose the distinction between t and l n , and 
so will at first obtain only the quotient curve ~K(N)/(wn)- We shall determine 
its double cover X(/V) by locating the branch points, which are the CM points 
on X(N)/ (wn) for which A is isomorphic to the product of two elliptic curves 
with CM by the quadratic imaginary order of discriminant — N or —47V; the 
arithmetic behavior of other CM points will then pin down the cover, including 
the right quadratic twist over Q. 

An abelian surface with QM by O has at least one principal polarization, 
and the number of principal polarizations of a generic surface with QM by O 
was computed in [Rotl, Theorem 1.4 and §6] in terms of the class number of 
Q(V _ AT). Each of these yields a map from X(N)/(wn) to A%, the moduli three- 
fold of principally polarized abelian surfaces. This map is either generically 1 : 1 
or generically 2:1, and in the 2 : 1 case it factors through an involution Wd = Wd> 
on X(7V)/(u>7v) where d, d' > 1 are integers such that N — dd' and 

(See the last paragraph of [Rot 21 §4], which also notes that a 2 : 1 map occurs 
for N = 6 and N = 10, each of which has a unique choice of polarization. In 
the other cases N = 14, 57, 206 that we study in this paper, only 1 : 1 maps 
arise, because the criterion (fTJ) is not satisfied.) We aim to determine at least one 
of the maps X(A r ) / (wn)—>A2 in terms of the Clebsch-Igusa coordinates on A2, 
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and thus to find the moduli of the generic abelian surface with endomorphisms 
by 00 

K3 surfaces, elliptic K3 surfaces, and the Dolgachev Kumar cor- 
respondence. Let F be a field of characteristic zero. Recall that a K3 surface 
over F is a smooth, complete, simply connected algebraic surface S/F with 
trivial canonical class. The Neron-Severi group NS(S') = NSp(S) is the group 
of divisors on S defined over the algebraic closure F, modulo algebraic equiv- 
alence. For a K3 surface this is a free abelian group whose rank, the Picard 
number p — p(S), is in {1, 2, 3, ... , 20}. The intersection pairing gives NS(S) the 
structure of an integral lattice; by the index theorem for surfaces, this lattice has 
signature (1, p — 1), and for a K3 surface the lattice is even: v ■ v = mod 2 for 
all v G NS(5). Over C, the cycle class map embeds NS(-S) into the "K3 lattice" 
H 2 (S,Z) = II 3 ,i 9 = U 3 8 E s (-lf, where U = ii M is the "hyperbolic plane" 
(the indefinite rank-2 lattice with Gram matrix (?j)), and Eg(— 1) is the E$ 
root lattice made negative-definite by multiplying the inner product by —1. The 
Torelli theorem of Piateckii-Shapiro and Safarevic [PSSj describes the moduli of 
K3 surfaces, at least over C: the embedding of NS(S') into 113,19 is primitive, 
that is, realizes NS(S') as the intersection of 113,19 with a Q-vector subspace of 
Il3,i9 ® Q; for every such lattice L of signature (l,p — 1), there is a nonempty 
(coarse) moduli space of pairs (S, t), where i : L — > NS(5) is a primitive embed- 
ding consistent with the intersection pairing; and each component of the moduli 
space has dimension 20 — p. Moreover, for p = 20, 19, 18, 17 these moduli spaces 
repeat some more familiar ones: isogenous pairs of CM elliptic curves for p = 20, 
elliptic and Shimura modular curves for p = 19, moduli of abelian surfaces with 
real multiplication or isogenous to products of two elliptic curves for p — 18, 
and moduli of abelian surfaces for certain cases of p = 17. Note the consequence 
that an algebraic family of K3 surfaces in characteristic zero with p > 19 whose 
members are not all F-isomorphic must have p = 19 generically, else there would 
be a positive-dimensional family of K3 surfaces with p > 20. 

An elliptic K3 surface S/F is a K3 surface together with a rational map 
t : S^P 1 , defined over F, whose generic fiber is an elliptic curve. The classes of 
the zero-section sq and fiber / in NS(S') then satisfy sq ■ sq = —2, sq ■ f = 1, and 
/ • / = 0, and thus generate a copy of U in NS(S') defined over F. Conversely, 
any copy of U in NS(5) defined over F yields a model of S as an elliptic surface: 
one of the standard isotropic generators or its negative is effective, and has 2 
independent sections, whose ratio gives the desired map to P 1 . We often use 
this construction to transform one elliptic model of S to another that would 
be harder to compute directly. (Warning: in general one might have to subtract 
some base locus from the effective generator to recover the fiber class /.) 

Since disc([/) = —1 is invertible, we have NS(5) = (so,f) © (so:f) ± , with 
the orthogonal complement (so,/) having signature (0, p — 2); we thus write 
( s 0i/) = N css (— 1) for some positive-definite even lattice N ess , the "essential 

2 Alas we cannot say simply "find the generic abelian surface with endomorphisms 
by O" , even up to quadratic twist, because there are abelian surfaces with rational 
moduli but no model over Q. 
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lattice" of the elliptic K3 surface. A vector v £ iV css of norm 2, corresponding to 
v G (so, /) X with v-v = —2, is called a "root" of N css ; let R C N css be the sublat- 
tice generated by the roots. This root sublattice decomposes uniquely as a direct 
sum of simple root lattices A n (n > 1), D n (n > 4), or E n (6 < n < 8). These 
simple factors biject with reducible fibers, each factor being the sublattice of N ess 
generated by the components of its reducible fiber that do not meet sq. The graph 
whose vertices are these components, and whose edges are their intersections, is 
then the A n , D n , or E n root diagram; if the identity component and its intersec- 
tion^) are included in the graph then the extended root diagram A n , D n , or E n 
results. The quotient group N ess /R is isomorphic with the Mordell-Weil group 
of the surface over F(t); the isomorphism takes a point P to the projection of the 
corresponding section s P to (s , f) 1 ' , and the quadratic form on the Mordell-Weil 
group induced from the pairing on N csa is the canonical height. Thus the Mordcll 
Weil regulator is r 2 disc(A^ 0SS ) /disc(i?) = r 2 |disc(NS(S))| / disc(E), where r is 
the size of the torsion subgroup of the Mordell-Weil group. 

An elliptic surface has Weierstrass equation Y 2 = X 3 + A(t)X + B{t) for 
polynomials A, B of degrees at most 8, 12 with no common factor of multiplicity 
at least 4 and 6 respectively, and such that either deg(A) > 4 or deg(-B) > 6 (i.e., 
such that the condition on common factors holds also at t — oo when A, B are 
considered as bivariate homogeneous polynomials of degrees 8, 12). The reducible 
fibers then occur at multiple roots of the discriminant A = — 16(4A 3 + TIB 2 ) 
where B does not vanish to order exactly 1 (and at t = oo if deg A < 22 and 
degi? ^ 11). To obtain a smooth model for S we may start from the surface 
Y 2 = X 3 + A(t)X + B(t) in the P 2 bundle P(O(0) 0(2) © 0(3)) over P 1 with 
coordinates (1 : X : Y), and resolve the reducible fibers, as exhibited in Tate's 
algorithm [Taj , which also gives the corresponding Kodaira types and simple root 
lattices. This information can then be used to calculate the canonical height on 
the Mordell-Weil group, as in [Si] . 

The Kummer surface Km(A) of an abelian surface A is obtained by blowing 
up the 16 = 2 4 double points of A/{±1}, and is a K3 surface with Picard 
number p{Km(A)) — p(A) + 16 > 17. In general NS(Km(A)) need not consist of 
divisors defined over F, even when NS(A) does, because each 2-torsion point of A 
yields a double point of A/{±1} whose blow-up contributes to NS(Km(A)), and 
typically Gal(F/F) acts nontrivially on A[2]. But when A is principally polarized 
Dolgachev |Doj constructs another K3 surface Sa/F, related with Km(A) by 
degree-2 maps defined over F, together with a rank-17 sublattice of NS(S'a) that 
is isomorphic with U © E-j © Es and consists of divisor classes defined over F. It 
is these surfaces that we parametrize to get at the Shimura curves X(A^). 

If A has QM then p(A) > 3, with equality for non-CM surfaces, so p(Sa) = 
p(Km(A)) > 19. When A has endomorphisms by O, we obtain a sublattice 
Ln C NS(SU) of signature (1, 18) and discriminant 2A^. This even lattice Ln is 
characterized by its signature and discriminant together with the following con- 
dition: for each odd p\N the dual lattice L* N contains a vector of norm c/p for 
some c G Z such that x P (c) = ~\ p {—2N/p), where \p ls the Legendre symbol 
{■/p) \ equivalently, N* ss contains a vector of norm c/p with x p (c) = — x P (+2A r /p). 
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There is a corresponding local condition at 2, but it holds automatically once 
the conditions at all odd p\N are satisfied; likewise when N is odd it is enough to 
check all but one p\N. The Shimura curve ~X(N)/(wn) parametrizes pairs («S, t) 
where S is a K3 surface with p(S) > 19 and I is an embedding Lat c — > NS(5). If 
= 20 then (5, t) corresponds to a CM point on X.(N)/(wn) whose discrim- 
inant equals disc(NS(S)). The CM points of discriminant — N or — AN are the 
branch points of the double cover X(JV) of X(AT) / (wn). The arithmetic of other 
CM points then determines the cover; for instance, if X(N)/ (wn) is rational, we 
know X(JV) up to quadratic twist, and then a rational CM point of discriminant 
D 7^ —TV, — AN lifts to a pair conjugate over Q(\/— D). 

The correspondence between and iSa was made explicit by Kumar |Ku|, 
Theorem 5.2]. Let A be the Jacobian of a genus-2 curve C, and let J2, Ia, Iq, ho 
be the Clebsch-Igusa invariants of C. (If a principally polarized abelian surface A 
is not a Jacobian then it is the product of two elliptic curves, and thus cannot 
have QM unless it is a CM surface.) We give an elliptic model of Sa with 
iV ess = R = E7 © Eg, using a coordinate ionP 1 that puts the E-j and Eg, fibers 
at t = and t = oo. Any such surface has the formula 

Y 2 = X 3 + (at 4 + a't 3 )X + (b"t 7 + bt 6 + b't 5 ) (2) 

for some a, a', 6, 6', b" with a', 6" 7^ 0. (There are five parameters, but the moduli 
space has dimension only 5 — 2 = 3 as expected, because multiplying t by a 
nonzero scalar yields an isomorphic surface, and multiplying a, a' by A 2 and 6, b' 
by A 3 for some A ^ yields a quadratic twist with the same moduli.) Kumar 
shows that setting 

(a, a', b, b\ b") = (-I4/I2, -1, {hh - 3/ 6 )/108, / 2 /24, J 10 /4) (3) 

in ([2|) yields the surface Sjm)- Starting from any surface (J2j) we may scale 
(t, X, Y) to (—a't, a' 2 X, a' 3 Y) and divide through by a' 6 to obtain an equation 
of the same form with a' = — 1; doing this and solving §S§ for the Clebsch-Igusa 
invariants Ii, we find 

(I 2 ,h,I 6 Jw) = {-24b' /a', -12a, 96ab'/a' - 36o, -Aa'b"). (4) 

If A has QM by 0, but is not CM, then the elliptic surface © has a Mordell- 
Weil group of rank 2 and regulator N, with each choice of polarization of A 
corresponding to a different Mordell-Weil lattice. The polarizations for which the 
map X(N)/(wn)— >Ai factors through some Wd are those for which the lattice 
has an involution other than — 1. When this happens, two points on X(N)/(wn) 
related by Wd yield the same surface © but a different choice of Mordell-Weil 
generators. For example, when N = 6 and N = 10 these lattices have Gram 
matrices | (® g) and | (® °) respectively. 

Some computational tricks. Often we need elliptic surfaces with an A n 
fiber for moderately large n, that is, for which AA 3 + 27 B 2 vanishes to moderately 
large order n + 1 at some t = to at which neither A nor _B vanishes. Thus we 



8 Noam D. Elkics 



have approximately (A, B) — (—3a 2 , 2a 3 ) near t = to- Usually one lets a be a 
polynomial that locally approximates (—A/3) 1 / 2 at t = to, and writes 

(A, B) = (-3(a 2 + 26), 2(a 3 + 3a6) + c) (5) 

for some b, c of valuations v(b) — v, v(c) = 2v at to- Then v(A) > 2v always, 
and v(A) > Zv if and only if v(36 2 — ac) > Sv; also if fj, < v then v(A) = 2v + p, 
if and only if v(3b 2 — ac) = 2v + fi. See |Haj : this was also the starting point of 
our analysis in [E2j . For our purposes it is more convenient to allow extended 
Weierstrass form and write the surface as 

Y 2 = X 3 + a(t)X 2 + 2b(t)X + c(t) (6) 

with polynomials a, b, c of degrees at most 4, 8, 12 such that (v(b), v(c)) — [v, 2v). 
Translating X by —a/3 shows that this is equivalent to with a, 6 divided 
by 3 (so /.i = v(b 2 — ac) in ©)• But © tends to produce simpler formulas, both 
for the surface itself and for the components of the fiber, which are rational if 
and only if a is a square. For instance, the Shioda-Hall surface with an Aig fiber 
[Sh Ha] can be written simply as 

Y 2 = X 3 + (t 4 + M 3 + 6t 2 + 7t + 4)X 2 - 2(< 3 + 2t 2 + 3t + 2)X + (t 2 + t + 1) 

with the A±g fiber at infinity, and this is the quadratic twist that makes all of 
NS(S') defined over Q. The same applies to D n , when A' := A/t 2 and B' := B/t 3 
are polynomials such that 4A' 3 + 27B' 2 has valuation n — 4. See for instance 
(|19[) below. When we want singular fibers at several t values we use an extended 
Weierstrass form (|6]) for which (v(b),v(c)) — (y,2v) holds (possibly with differ- 
ent v) at each of these t. 

Having parametrized our elliptic surface S with Ltv c ~ > NS(S'), we seek spe- 
cializations of rank 20 to locate CM points. In all but finitely many cases S has 
an extra Mordell-Weil generator. In the exceptional cases, either some of the re- 
ducible fibers merge, or one of those fibers becomes more singular, or there is an 
extra A\ fiber. Such CM points are easy to locate, though some mergers require 
renormalization to obtain a smooth model and find the CM discriminant D, as 
we shall see. When there is an extra Mordell-Weil generator, its height is at 
least \D\/2N, but usually not much larger. (Equality holds if and only if the ex- 
tra generator is orthogonal to the generic Mordell-Weil lattice; in particular this 
happens if S has generic Mordell-Weil rank zero.) The larger the height of the 
extra generator, the harder it typically is to find the surface. This has the curious 
consequence that while the difficulty of parametrizing S increases with N, the 
CM points actually become easier to find. In some cases we cannot solve for the 
coefficients directly. We thus adapt the methods of |E3j . exhaustively searching 
for a solution modulo a small prime p and then lifting it to a p-adic solution 
to enough accuracy to recognize the underlying rational numbers. We choose 
the smallest p such that Xp(~ D) = +1, so that reduction mod p does not raise 
the Picard number, and we can save a factor of p in the exhaustive search by 
first counting points mod p on each candidate S to identify the one with the 
correct CM. 
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For large N we use the following variation of the p-adic lifting method to 
find the Shimura curve X(iV)/ (wn) and the surfaces S parametrized by it. First 
choose some indefinite primitive sublattice L' C Ln and parametrize all S with 
NS(S') 3 V . Search in that family modulo a small prime p to find a surface 
So with the desired L^- Let /i,/2 be simple rational functions on the (S,L') 
moduli space. We hope that the degrees, call them di, of the restriction of fi 
to X(iV)/ (wn) are positive but small; that /i is locally 1 : 1 on the point of 
X(N)/(wn) parametrizing So] and that the map (fi,fz) ■ X(N)/(wn)~ >A 2 is 
generically 1 : 1 to its image in the affine plane. For various small lifts /i of 
fi(So) to Q, lift So to a surface S/Q p with fx(S) — fi, compute /^(S 1 ) to high 
p-adic precision, and use lattice reduction to recognize f2(S) as the solution of 
a polynomial equation -FX/2) — of degree (at most) d\. Discard the few cases 
where the degree is not maximal, and solve simultaneous linear equations to 
guess the coefficients of F as polynomials of degree at most d 2 in f\. At this 
point we have a birational model F(fi, f%) = for X(N)/(wn}- Then recover a 
smooth model of the curve (using Magma if necessary) , recognize the remaining 
coefficients of S as rational functions by solving a few more linear equations, and 
verify that the surface has the desired embedding NS(S). 

3 N = 6: The first Shimura curve 

The K3 surfaces. We take N css — R = A 2 © D-j © Eg, which has discriminant 
3 • 4 • 1 = 12 = 2N, and the correct behavior at 3 because A 2 contains vectors of 
norm 2/3 with x 3 (2) = — Xa(2 ■ 6/3)[= — 1]. We choose the rational coordinate t 
on P 1 such that the A 2 , D-j, and E$ fibers are at t = 1, 0, and oo respectively. 
If we relax the condition at t — 1 by asking only that the discriminant vanish to 
order at least 2 rather than 3 then the general such surface can be written as 

Y 2 = X 3 + (ao + a^tX 2 + 2a bt 3 {t - l)X + a Q b 2 t 5 (t - l) 2 (7) 

for some ao,ai,b, with a\b ^ lest the surface be too singular at t = 0. The 
discriminant is then t 9 (t— l) 2 Ai(t) with A\ a cubic polynomial such that A(l) = 
— 64aoai(ao + ai) 2 6 2 . Thus Ai(l) = if and only if a\ = or ao + a\ = 0. In 
the latter case the surface has additive reduction at t — 1. Hence we must 
have a\ — 0. The non-identity components of the resulting A 2 fiber at t = 1 
then have X = 0(t — 1); we calculate that X = x%(t — 1) + 0((t — l) 2 ) makes 
Y 2 = (xi + b) 2 ao(t — l) 2 + 0(t — l) 3 . Therefore these components are rational 
if and only if ao is a square. We can then replace {X, Y, b) by (a X, a^ 2 Y, a b) 
in (J7]) to obtain the formula 

Y 2 = X 3 + tX 2 + 2bt 3 (t - l)X + b 2 t 5 (t - l) 2 (8) 

for the general elliptic K3 surface with N css = R = A-iD^Eg and rational A 2 
components. The two components of the D-j fiber farthest from the identity com- 
ponent then have X = bt 2 + 0(t 3 ), so Y 2 = b 3 t 6 + 0(t 7 ); thus these components 
are both rational as well if and only & is a square, say b — r 2 . Then b and r 
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are rational coordinates on the Shimura curves X(6)/(i02> ^3) and ~X.(6)/(we) 
respectively, with the involution u>2 = u>3 on ~X.(6)/(we) taking r to — r. 

The elliptic surface © has discriminant A = 16bH 9 (t - l) 3 (276(i 2 - f) - 4). 
Thus the formula (JSJ fails at 6 — 0, and also of course at b = 00. Near each of 
these two points we change variables to obtain a formula that extends smoothly 
to b = or b = 00 as well. These formulas require extracting respectively a 
fourth and third root of (3, presumably because 6 = and b — 00 are elliptic 
points of the Shimura curve. For small b, we take b — (3 4 and replace (t, X, Y) 
by (t/(3 2 ,X/[3 2 ,Y/(3 3 ) to obtain 

Y 2 = X 3 + tX 2 + 2t 3 (t - (3 2 )X + t 5 (t - (3 2 ) 2 , (9) 

with the A2 fiber at t — (3 2 rather than t = 1. When f3 = 0, this fiber merges 
with the Dj fiber at t = to form a -D10 fiber, but we still have a K3 surface, 
namely Y 2 — X 3 + tX 2 + 2t 4 X + t\ with L = R = D w © Eg. This is the CM 
point of discriminant —4. For large b, we write b — 1//3 3 and replace (X,Y) by 
(X//3 2 ,F//3 3 ) to obtain 

Y 2 = X 3 + (3 2 tX 2 + 2/3t 3 (t - 1)X + t 5 (t~ l) 2 ; (10) 

then taking fi^O yields the surface Y 2 — X 3 + t 5 (t - l) 2 with A^ css = R = 
A2 © E$ Eg: the t = fiber changes from D7 to £«, and the t = 1 fiber 
becomes additive but still contributes A2 to i? (Kodaira type IV rather than I3). 
This is the CM point of discriminant —3. 

Two more CM points. The factor 27b(t 2 — t) — 4 of A is a quadratic 
polynomial in t of discriminant 27b(27b + 16). Hence at b = —16/27 we have 
iV ess = R = Ai(BA2®Dt(BEs, and we have located the CM point of discriminant 
—24. Three points fix a rational coordinate on P 1 , so we can compare with the 
coordinate used in |E1( Table 1], which puts the CM points of discriminant —3, 
—4, and —24 at 00, 1, and respectively; thus that coordinate is 1 + 276/16. This 
also confirms that X(6) is obtained by extracting a square root of — (27r 2 + 16). 

We next locate a CM point of discriminant —19 by finding b for which the 
surface (|8]) has a section s p of canonical height 19/12. This is the smallest 
possible canonical height for a surface with R = A2 © -D7 © E$, because the 
naive height is at least 4 and the height corrections at the A2 and D7 fibers can 
reduce it by at most 2/3 and 7/4 respectively, reaching 4 — 2/3 — 7/4 = 19/12. 
Let (X(t), Y{t)) be the coordinates of a point P of height 19/12. Then X(t) and 
Y(t) are polynomials of degree at most 4 and 6 respectively (else s P intersects So 
and the naive height exceeds 4), and X vanishes at t = 1 (so s P passes through 
a non-identity component of the A2 fiber) and has the form bt 2 + 0(t 3 ) at t = 
(so s P meets one of the components of the D7 fiber farthest from the identity 
component). That is, X — b(t 2 — t 3 )(l + tit) for some t\. Substituting this into 
© and dividing by the known square factor (t 4 — t 3 ) 2 yields b 3 times 

- tit 4 + (t 3 - 3t 2 )t 3 + 3(t 2 - h)t 2 + ((3ti - 1) + b-H 2 )t + 1, (11) 

so we seek 6, t\ such that the quartic (fTTj) is a square. We expand its square root 
in a Taylor expansion about t — and set the t 3 and t 4 coefficients equal to 
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zero. This gives a pair of polynomial equations in 6 and t\, which we solve by 
taking a resultant with respect to t\. Eliminating a spurious multiple solution 
at b = 0, we finally obtain (b,tx) — (81/64,-9), and confirm that this makes 
(fTTj) a square, namely (27t 2 — 18t— l) 2 . Therefore 81/64 is the ^-coordinate of a 
CM point of discriminant —19. Then 1 + 276/16 = 3211/2 10 , same as the value 
obtained in IE1I. 



Clebsch— Igusa coordinates. The next diagram shows the graph whose 
vertices are the zero-section (circled) and components of reducible fibers of an 
elliptic K3 surface S with N css = A2®Di®E^, and whose edges are intersections 
between pairs of these rational curves on the surface. Eight of the vertices form 
an extended root diagram of type EV, and are marked with their multiplicities in 
a reducible fiber of type Ej of an alternative elliptic model for S. We may take 
cither of the unmarked vertices of the D-? subgraph as the zero-section. Then 
the essential lattice of the new model includes an Eg root diagram as well as the 
forced £7. We can thus apply Kumar's formulas to this model once we compute 
its coefficients. 



D 7 


1 2 3 


2 

4 3 






Figure 1: An £7 divisor supported on the zero-section and 
fiber components of an AiD^Es surface 

The sections of the £7 divisor are generated by 1 and u := Xj (t 4 — t 3 ) + b/t. 
Thus u : S^V 1 gives the new elliptic fibration. Taking X = (t 3 ~t 2 )(tu-b) in (JSJ 
and dividing by (i 4 - 1 3 ) 2 yields Y 2 = Q(t) for some quartic Q. Using standard 
formulas for the Jacobian of such a curve, and bringing the resulting surface 
into Weierstrass form, we obtain a formula ^ with (a, a' ,b,b' ,b") replaced by 
(-36, 1, -2b 2 , -(b + 1), -b 3 ). As expected this surface has Mordell-Weil rank 2 
with generators of height 5/2, namely 

(r 6 t 4 + 2(r i + r 3 )t 3 + (r 2 + l)t 2 , r 9 t 6 + 3(r 7 + r 6 )t 5 +3(r 5 +r 4 + r 3 )t A + {r 3 + l)t 3 ) 



and the image of this section under r^ — r (recall that 6 = r 2 ). The formula ([4]) 
yields the Clebsch-Igusa coordinates 



(I 2 ,h,I 6 ,I w ) = ((246+ 1), 366, 726(56 + 4), 46 3 ). 



(12) 
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4 N = 14: The CM point of discriminant —67 

The K3 surfaces. Here we take 7V CSS = R = A3 A§ Eg, which has discrim- 
inant 4-7-1 = 28 = 2N, and the correct behavior at 7 because A§ contains 
vectors of norm 6/7 with Xri®) = ~Xt(2 ■ 14/7)[= — 1]. We put the A3, A 6 , and 
E 8 fibers at t = 1,0, and 00 respectively. We then seek an extended Weier- 
strass form © with a,b,c of degrees 2,4,7 such that t 3 — t 2 \b, (t 3 — t 2 ) 2 |c, 
and (t 3 — t 2 ) 6 \b 2 — ac. This gives at least A3, A^, Eg. It is then easy to impose 
the extra condition t 7 \A, and we obtain a — A((s + l)i 2 + (3s 2 + 2s)t + s 3 ), 
b = X 2 (s + l)((4s + 2)t + 2s 2 )(t 3 - t 2 ), c = A 3 (s + l) 2 (t + s){t 3 - t 2 ) 2 for some 
s, A. The twist A must be chosen so that a(0) and a(l) are both squares; this is 
possible if and only if s 2 + s is a square, so s = r 2 /(2r + 1) for some r. Thus r and 
s are rational coordinates on X(14)/(u>i4) and X(14)/(u>2, wr) respectively, with 
the involution W2 — W7 on X(14)/(wi4) taking r to ~r/(2r + 1). The formula 
in terms of r is cleaner if we let the A3 fiber move from t = 1; putting it at 
t = 2r + 1 yields 

a = ((r + 1) 2 )< 2 + (3r 4 + 4r 3 + 2r 2 )t + r 6 ), 

b = 2(r + l) 2 ((2r 2 + 2r + l)t + r 4 )(t - (2r + l))t 2 , (13) 
c= (r + l) 4 (t-(2r + l)) 2 (t + r 2 )t 4 . 

Easy CM points. At r = 0, the A§ fiber becomes E-j, so we have a CM 
point with D = —8; at r = —1/2, the A3 and ^6 fibers merge to A\q, giving a 
CM point with D = —11. These have s = 0, s = 00 respectively. There is an 
extra Ai fiber when lis 2 + 3s + 8 = 0; the roots of this irreducible quadratic 
give the CM points with D = —56 (and their lifts to X(14)/ (11)14) are the branch 
points of the double cover X(14)). In [Elj we gave a rational coordinate t on 
X(14) / (w2,wr) for which the CM points of discriminants —8, —11, and —56 had 
t = 0, t = —1, and 16i 2 + 13< + 8 = respectively. Therefore that t is our 

-8/(8 + 1). 

A harder CM point. At the CM point of discriminant —67 our surface has 
a section of height 67/28 = 4 - (3/4) - (6/7). Thus Y 2 = X 3 + aX 2 + bX + c 
has a solution in polynomials X, Y of degrees 4, 6 with X(0) — X(2h + 1) = 
and Y having valuation exactly 1 at t = and t = 2h + 1. An exhaustive search 
mod 17 quickly finds an example, whose lift to Q17 then yields r = —35/44 with 

X = -p-^ t [22t + 13) (527076£ 2 + 760364i + 275625). 

Thus s = -1225/1144, and -s/(s + 1) confirms the entry -1225/81 in the 
\D\ = 67 row of [ED Table 5]. 

5 N — 57: The first curve ~K(N)/wn of positive genus 

The K3 surfaces. We cannot have iV css = R here because there is no root 
lattice of rank 17 and discriminant 6-19. Instead we take for R the rank-16 lattice 
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A$ © An of discriminant 6 ■ 12 = 72, and require an infinite cyclic Mordell-Weil 
group N oss /R with a generator corresponding to a section that meets the A$ and 
An fibers in non-identity components farthest from the A§ identity and nearest 
the An identity respectively, and does not meet the zero-section (i.e., for which 
X is a polynomial of degree at most 4 in t) . Such a point has canonical height 



3 • 3 



1 • 11 
12 



19 
12 



2N 
disc R 



(14) 



Thus disc(-/V ess ) has the desired discriminant 27V. We may check the local con- 
ditions by noting that A§ contains a vector of norm 4/3 that remains in AT* SS , 
and x 3 (4) = -X3C 2 • 57/3) [= +1]. We put the A 5 fiber at t = and the An 
fiber at t = 00. We eventually obtain the following parametrization in terms of 
a coordinate r on the rational curve X(57)/ (1V3, wig): let 



p(r)=4(r-l)(r 2 -2) + l, 

d= (r 2 - l) 2 (9t+ (2r- l)p(r)), 

c = 9t 2 - (2r - l)(8r 2 + 4r - 22)< + (2r - l) 2 p(0> 

& = (t- (r 2 - 2r))c + d, 

2 o ™ \ \ 2 „ 1 o/'-/- /^.2 o,v,\ \ ^7 1 /™2 



(15) 



a=(t- (r 2 - 2r)) 2 c + 2{t - (r 2 - 2r))d + (r 2 - l) 4 ((4r + A)t + p(r)); 
Then the surface is 



Y 2 = X A + aX 2 + 8(r - l) 4 (r + l) 5 6t 2 X + 16(r - l) 8 (r + l) 10 ci 
with a section of height 19/12 at 

4(r- l) 4 (r + l) 5 (2r- l)t 2 



4+3 



X = -- 



(r 2 



1) 



4(r-2)(r + l) 4 i 
r 2 — r + 1 



(16) 



(17) 



The components of the An fiber are rational because the leading coefficient 
of a is 9, a square; the constant coefficient is (r 2 — r + l) 4 p(r), so ~K(57)/(w5?) 
is obtained by extracting a square root of p(r). This gives the elliptic curve 
with coefficients [ai, a 2 , 03, a 4 , ag] = [0,-1,1,-2,2], whose conductor is 57 as 
expected (see e.g. Cremona's tables [Crj where this curve appears as 57-Al(E)). 

This curve has rank 1, with generator P = (2, 1). The point at infinity is the 
CM point of discriminant —19; this may be seen by substituting l/s for r and 
(t/s 3 ,X/s 12 ,Y/s ls ) for (t,X,Y), then letting s^O to obtain the surface 



Y 2 = X J 



(9t 4 - 16t 3 + U)X 2 + {72t 5 - l28t 4 )X + (144t 6 - 256i 5 



(18) 



with a Dq fiber at t = rather than an A5. Then we still have a section (X, Y) = 
(4i 3 -8i 2 , (3i-5)(t 4 -t 3 )) of height 19/12, but there is a 2-torsion point (X,Y) = 
(-4t,0) so disc(7V 0SS ) = -disc(NS(S')) drops to 4 • 12 • (19/12)/2 2 = 19. The 
remaining rational CM points on X(57)/(u?57) come in six pairs ±nP: 



n 


1 


2 


3 


4 


5 


8 


r 


2 


1 


-1 





5/4 


13/9 


-D 


7 


4 


16 


28 


43 


163 
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The last three of these have extra sections X — — At, X — 0, and 

X = -28 • ll 3 (0 2 /3 6 ) + (415454</3 18 )) 

respectively. At r = 2, the v4 n fiber becomes an A 12 and our generic Mordell- 
Weil generator becomes divisible by 3; the new generator (— 972i, 26244i 2 ) has 
height 4 - (5/6) - (40/13) = 7/78, so disc iV 0SS = 6 • 13 • (7/78) = 7. At r = 1, the 
A5 and An fibers together with the section all merge to form a Dig fiber: let 
r = 1 + s and change (t, X) to (st — 1, — 8s 3 X), divide by (—2s) 9 , and let s— >0 
to obtain the second Shioda-Hall surface 

X 3 + (t 3 + 8t)X 2 - (32< 2 + 12S)X + 256t (19) 

with a Dis fiber at t = 00 [ShIHaj . At t = — 1, the reducible fibers again merge, 
this time forming an An while the Mordell-Weil generator's height drops to 
4 - (4 • 14/18) = 8/9, whence disc(7V oss ) = 16. 

We find four more rational CM values of r that do not lift to rational points 
on X(57)/(ui 57 ), namely r = 5, 1/2, 17/16, —7/4, for discriminants —123, —24, 
—267, and —627 = —11-57 respectively. The first of these again has an A12 fiber, 
this time with the section of height 4 — (9/6) — (12/13) = 41/26; the second has 
a rational section at X = 0; in the remaining two cases we find the extra section 
by p-adic search: 

ll3o2 

X = - 22lQi ^ 2 (7840t 2 - 2037t + 3267) (20) 

for r = 17/16, and 

x = 3W g (t) 

2 12 (81920i 3 + 9216i 2 + 23868t + 39339) 2 V ' 

for r — —7/4, where q(t) is the quintic 

419430400i 5 + 2846883840£ 4 + 17148174336* 3 

+ 78784560576i 2 + 175272616341* - 12882888. 

Using I we can show that there are no further rational CM values. 



6 N — 206: The last curve ~K(N)/w N of genus zero 

Summary of results. Again we take AT ess of rank 16 and an infinite cyclic 
Mordell-Weil group, here R = A 2 © At © A w with a Mordell-Weil generator of 
height 412/165 = 6 - (1 • 2/3) - (2 • 3)/5 - (2 • 9)/ll. With the reducible fibers 
placed at 1,0, 00 as usual, the choice of R means A = t 5 (t — l) 3 Z\i with Ai of 
degree 24— (3 + 5 + 11) = 5 and ^i(0), ^i(l) ^ 0; the generator must then have 
X(t) = X\(t)/(t — to) 2 for some sextic Xi and some to 7^ 0, 1, with the corre- 
sponding section passing through a non-identity component of the A 2 fiber and 
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components at distance 2 from the identity of A4 and A w . We eventually suc- 
ceed in parametrizing such surfaces, finding a rational coordinate on the modular 
curve X(206)/(u>206)- These elliptic models do not readily exhibit the involution 
W2 = W103 on this curve, so we recover this involution from the fact that it 
must permute the branch points of the double cover X(206) of X(206)/(u>2u6)- 
We locate these branch points as simple zeros of the discriminant of A\. As 
expected, there are 20 (this is the class number of Q(\/— 206)), forming a single 
Galois orbit. We find a unique involution of the projective line X(206)/(ii;206) 
that permutes these zeros. This involution has two fixed points, so we switch 
to a rational coordinate r on X(206)/(u;206) that makes the involution r^—r. 
Then ro := r 2 is a rational coordinate on X(206)/ (1V2, W103), and the 20 branch 
points are the roots of Pin(r 2 ) where Pio is the degree-10 polynomial 

Pw(r ) = 8ro° - 13rg - 42r* - 331^ - 220r[j + 733^ (22) 
+ 6646r^ + 19883rj] + 28840^ + 18224r + 4096. 

As a further check on the computation, Pio has dihedral Galois group, discrim- 
inant — 2 138 1 03 7 , and field discriminant — 2 12 103 5 , while Pio(r 2 ) has discrimi- 
nant 2 311 103 14 and field discriminant 2 27 103 10 . We find that r = 0,±l,±2,oo 
give CM points of discriminants D = —4, —19, —163, —8 respectively; evaluating 
Pw(t 2 ) at any of these points gives — D times a square, showing that the Shimura 
curve X(206) has the equation s 2 = —Pio(r 2 ) over Q. The curves X(206)/(w2), 
X(206)/(wio3) are then the double covers s 2 , = —Pio(ro), s 2 = — roPio(^o) of 
the ro-line X(206)/(w2, W103) ( m that order, because W103 cannot fix a CM point 
of discriminant —4 or —8). 
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